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Abstract. By a combination of asymptotic ODE estimates and nu- 
merical Evans function calculations, we establish stability of vis- 
cous shock solutions of the isentropic compressible Navier-Stokes 
equations with 7-law pressure (i) in the limit as Mach number M 
goes to infinity, for any 7 > 1 (proved analytically), and (ii) for 
M > 2,500, 7 € [1,2.5] (demonstrated numerically). This builds on 
and completes earlier studies by Matsumura-Nishihara and Barker- 
Humpherys-Rudd-Zumbrun establishing stability for low and inter- 
mediate Mach numbers, respectively, indicating unconditional stabil- 
ity, independent of shock amplitude, of viscous shock waves for 7-law 
gas dynamics in the range 7 S [1, 2.5]. Other 7- values may be treated 
similarly, but have not been checked numerically. The main idea is to 
establish convergence of the Evans function in the high-Mach number 
limit to that of a pressureless, or "infinitely compressible" , gas with 
additional upstream boundary condition determined by a boundary- 
layer analysis. Recall that low-Mach number behavior is incompress- 
ible. 



1. Introduction 

The isentropic compressible Navier-Stokes equations in one spatial 
dimension expressed in Lagrangian coordinates take the form 



vt - 



= 




(1) 
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where v is specific volume, u is velocity, and p pressure. We assume 
an adiabatic pressure law 

p{v) = aov~^ (2) 

corresponding to a 7-law gas, for some constants ao > and 7 > 1. 
In the thermodynamical rarified gas approximation, 7 > 1 is the 
average over constituent particles of 7 = (n + 2)/n, where n is the 
number of internal degrees of freedom of an individual particle 
n = 3 (7 = 1.66...) for monatomic, n = 5 (7 = 1.4) for diatomic 
gas. For dense fluids, 7 is typically determined phenomenologically 
|19j . In general, 7 is usually taken within 1 < 7 < 3 in models of 
gas- or fluid-dynamical flow, whether phenomenological or derived 
by statistical mechanics [32 t [39 t [30j . 

It is well known that these equations support viscous shock waves, 
or asymptotically-constant traveling-wave solutions 

{v,u){x,t) = {v,u){x — st), lim {v,u){z) = {v,u)±, (3) 

2— > it 00 

in agreement with physically-observed phenomena. In nature, such 
waves are seen to be quite stable, even for large variations in pressure 
between v±. However, it is a long-standing mathematical question to 
what extent this is reflected in the continuum- mechanical model ([T]), 
that is, for which choice of parameters m±, 7) are solutions of ([3]) 
time-evolutionarily stable in the sense of PDE; see, for example, the 
discussions in [23113]. 

The first result on this problem was obtained by Matsumura and 
Nishihara in 1985 [34] using clever energy estimates on the integrals 
of perturbations in v and u, by which they established stability with 
respect to the restricted class of perturbations with "zero mass", i.e., 
perturbations whose integral is zero, for shocks with sufficiently small 
amplitude 

\p{v+) -p{v^)\ < C(t;_,7), 

with C — > 00 as 7 — > 1, but C « 00 for 7 7^ 1. There followed a 
number of works by Liu, Goodman, Szepessy-Xin, and others [28l[T5l 
H3|l29j toward the treatment of general, nonzero mass perturbations; 
see |47 l H5 t ,46j and references therein. A complete result of stability 
with respect to general n perturbations of small-amplitude 
shocks of system ([1]) was finally obtained in 2004 by Mascia and 
Zumbrun ^33j using pointwise semigroup techniques introduced by 
Zumbrun and Howard [471144] in the strictly parabolic case. 

The result of [33], together with the small-amplitude spectral sta- 
bility result of Humpherys and Zumbrun [22^1 generalizing that of 

^ The result of [21] is obtained by energy estimates combining the techniques of 
[34j with those of [141115] : a similar approach has been used in [30] to obtain small- 
amplitude zero-mass stability of Boltzmann shocks. See [371112] for an alternative 
approach based on asymptotic ODE methods. 
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[34j , in fact yields stability of small-amplitude shocks of general sym- 
metric hyperbolic-parabolic systems, largely settling the problem of 
small-amplitude shock stability for continuum mechanical systems. 

However, there remains the interesting question of large- amplitude 
stability. The main result in this direction, following a general strat- 
egy proposed in [l^; is a "refined Lyapunov theorem'tj established 
by Mascia and Zumbrun |32y46j for general symmetric hyperbolic- 
parabolic systems, stating that linearized and nonlinear n 
L^nH^ orbital stability (the standard notions of stability) are equiv- 
alent to spectral stability, or nonexistence of nonstable (nonnegative 
real part) eigenvalues of the linearized operator L about the wave, 
other than the single zero eigenvalue arising through translational 
invariance of the underlying equations. 

This reduces the problem of large-amplitude stability to the study 
of the associated eigenvalue equation {L — X)u = 0, a standard analyt- 
ically and numerically well-posed (boundary value) problem in ODE, 
which can be attacked by the large body of techniques developed for 
asymptotic, exact, and numerical study of ODE. In particular, there 
exist well-developed and efficient numerical algorithms to determine 
the number of unstable roots for any specific linearized operator L, 
independent of its origins in the PDE setting; see, e.g., [9l[T0l[m[8l[22] 
and references therein. In this sense, the problem of determining sta- 
bility of any single wave is satisfactorily resolved, or, for that matter, 
of any compact family of waves. To determine stability of a family 
of waves across an unbounded parameter regime, however, is another 
matter. It is this issue that we confront in attempting to decide the 
stability of general isentropic Navier-Stokes shocks. 

As pointed out in [i7ll2T| . zero-mass stability implies (and in a 
practical sense is roughly equivalent to) spectral stability. Thus, the 
original results of Matsumura and Nishida [3l] imply small-amplitude 
shock stability for general 7 and large-amplitude stability as 7 ^ 1. 
Recently, Barker, Humpherys, Rudd, and Zumbrun [2j have carried 
out a numerical Evans function study indicating stability on the large, 
but still bounded, parameter range 7 G [1,3], 1 < M < 3,000, 
where M is the Mach number associated with the shock. For dis- 
cussion of the Evans function, see Section 12.41 Recall that Mach 
number is an alternative measure of shock strength, with 1 corre- 
sponding to \p{v+) — p{v-)\ = and M — > 00 corresponding to 
\p{v+) — p{v-)\ — > 00; see Appendix A, [2J, or Section [2T] below. 
Mach 3, 000 is far beyond the hypersonic regime M ~ 10^ encoun- 
tered in current aerodynamics. However, the mathematical question 
of stability across arbitrary 7, M remains open up to now. 



"Refined" because the linearized operator L does not possess a spectral gap, 
hence e^' decays time-algebraically and not exponentially; see [471146] for further 
discussion. 
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In this paper, we resolve this question, using a combination of 
asymptotic ODE estimates and Evans function calculations to con- 
clude, first, stability of isentropic Navier-Stokes shocks in the large 
Mach number limit M oo for any 7 > 1, and, second, stability for 
ah M > 2, 500 for 7 G [1, 2.5] (for 7 E [1, 2], we obtain in fact stabil- 
ity for M > 500.) The first result is obtained analytically, the second 
by a systematic numerical study. Together with the numerical results 
of [2], this gives convincing numerical evidence of unconditional sta- 
bility for ^ € [1,2.5], independent of shock amplitude. As in [2], our 
numerical study is not a numerical proof, but contains the neces- 
sary ingredients for one; see discussion. Section [6j The restriction to 
7 G [1, 2.5] is an arbitrary one coming from the choice of parameters 
on which the numerical study [2] was carried out; stability for other 
7 can be easily checked as well. 

Our method of analysis is straightforward, though somewhat del- 
icate to carry out. Working with the rescaled and conveniently rear- 
ranged versions of the equations introduced in [2], we observe that the 
associated eigenvalue equations converge uniformly as Mach number 
goes to infinity on a "regular region" x < L, for any fixed L > 0, to 
a limiting system that is well-behaved (hence treatable by the stan- 
dard methods of |3i y371 [2]) in the sense that its coefficient matrix 
converges uniformly exponentially in x to limits at a; = ±00, but is 
underdetermined ai x = +00. 

On the complementary "singular region" x > L, the convergence 
is only pointwise due to a fast "inner structure" featuring rapid vari- 
ation of the converging coefficient matrices near x = +00, but the be- 
havior at X = -|-oo is of course determinate. Performing a boundary- 
layer analysis on the singular region and matching across x = L, we 
are able to show convergence of the Evans function of the original 
system as the Mach number goes to infinity to an Evans function of 
the limiting system with an appropriately imposed additional condi- 
tion at X = -|-oo, upstream of the shock. This reduces the question of 
stability in the high-Mach number limit to existence or nonexistence 
of zeroes of the limiting Evans function on 5ReA > 0, a question that 
can be resolved by routine numerical computation as in [2j, or by 
energy estimates as in Appendix [Bl 

The limiting system can be recognized as the eigenvalue equation 
associated with a pressure-less (7 = 0) gas, that is, the "infinitely- 
compressible" limit one might expect as the Mach number goes to 
infinity. Recall that behavior in the low Mach number limit is incom- 
pressible |251I271[T8] . However, the upstream boundary condition has 
to our knowledge no such simple interpretation. Indeed, to carry out 
the boundary-layer analysis by which we derive this condition is the 
main technical difficulty of the paper. 

Besides their independent interest, the results of this paper seem 
significant as prototypes for future analyses. Our calculations use 
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some properties specific to the structure of ([I]). In particular, we 
make use of surprisingly strong energy estimates carried out in [2] 
in confining unstable eigenvalues to a bounded set independent of 
shock strength (or Mach number). Also, we use the extremely simple 
structure of the eigenvalue equation to carry out the key analysis of 
the eigenvalue flow in the singular region near x = +00 essentially 
by hand. However, these appear to be conveniences rather than es- 
sential aspects of the analysis. It is our hope that the basic argument 
structure of this paper together with [2j can serve as a blueprint for 
the study of large- amplitude stability in more general situations. 

In particular, we expect that the analysis will carry over to the full 
(nonisentropic) equations of gas dynamics with ideal gas equation of 
state, which, formally, decouple in the high-Mach number limit into 
the equations of isentropic pressureless gas dynamics studied here, 
augmented with a separate temperature equation governed by simple 
diffusion/Fourier's law. 

2. Preliminaries 

We begin by recalling a number of preliminary steps carried out in [2]. 
Making the standard change of coordinates x ^ x — st, we consider 
instead stationary solutions {v,u){x,t) = {v,u){x) of 



where e is chosen so that < < = 1, our system takes the 
convenient form 



Vt — SVx — U; 



= 





Under the rescaling 



{x,t,v,u) — > {—esx,es'^t,v/e,—u/{es)), 



(5) 



= 





where a = uqE s 



2.1. Profile equation 



Steady shock profiles of ([6j) satisfy 

v' -u' = Q 
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subject to boundary conditions {v,u){±oo) = {v±,u±), or, simplify- 
ing, 



v' + {av-'')' = 



V 

Integrating from — oo to x, we get the profile equation 

v' = H{v,v+) ■.= v{v -l + aiv-"/ -1)), (7) 

where a is found by setting x = +oo, thus yielding the Rankine- 
Hugoniot condition 

^K-, IT- (8) 



V, 



Evidently, 0^7"^ in the weak shock limit ^ 1, while a ^ in 
the strong shock limit — > 0. The associated Mach number M may 
be computed as in Appendix A, as 

M = (7a)-i/2 (9) 

so that M ~ ^^^I'^v^^'^ +00 as f 4. — > and M ^ 1 as 11+ ^ 1; 
that is, the high-Mach number limit corresponds to the limit ^ 0. 



2.2. Eigenvalue equations 

Linearizing ([6|) about the profile {v,u), we obtain the eigenvalue prob- 
lem 



Xv + v' — u = 
Xu + u' ■ 



h{v) \' fu'\ (10) 
^ 1 _ I 1 



where 



/i(^)) = + a(7 - 1) + (a + l)z)^. (11) 

We seek nonstable eigenvalues A € {5Re(A) > 0} \ {0}, i.e., A for 
which ([To]) possess a nontrivial solution (f , u) decaying at plus and 
minus spatial infinity. As pointed out in |47ll21j . by divergence form 
of the equations, such solutions necessarily satisfy j'^^ v{x)dx = 
u{x)dx = 0, from which we may deduce that 

u{x) = I u{z)dz, v{x) = / v{z)dz 
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and their derivatives decay exponentially as a; ^ oo. Substituting and 
then integrating, we find that (u, v) satisfies the integrated eigenvalue 
equations (suppressing the tilde) 

Xv + v' -u = 0, (12a) 

, h(v) I u" , , s 

Xu + u - -:^v' = —. 12b) 

This new eigenvalue problem differs spectrally from (llOj) only at A = 
0, hence spectral stability of (fTOl) is equivalent to spectral stability of 
(jl2p . Moreover, since ()12p has no eigenvalue at A = 0, one can expect 
more uniform stability estimates for the integrated equations in the 
vicinity of A = [SllIllllZ!- 



2.3. Preliminary estimates 

The following estimates established in [2] indicate the suitability of 
the rescaling chosen in Section [2.11 For completeness, we prove these 
in Appendix [Xl 

Proposition 1 ([2]). For each 7 > 1, < u+ < 1, ([7]) has a unique 
(up to translation) monotone decreasing solution v decaying to its 
endstates with a uniform exponential rate, independent ofv-^., 7. In 
particular, for < w+ < ^ and v{0) := t>+ + 

\v{x) -v+\ < (^-^^e'^ x>0, (13a) 

\v{x) -v-\ < (^j^e^ x<0. (13b) 

Proposition 2 ([2]). Nonstable eigenvalues X of (jl2p . i.e., eigen- 
values with nonnegative real part, are confined for any < u+ < 1 to 
the region 

3f?e(A) + |9m(A)| < (^7 + ^)'- (14) 



2.4. Evans function formulation 

Following [2], we may express (fT2|) concisely as a first-order system 

W' = A{x,X)W, (15) 



Aix,X)=(oO 1 ], W=(l], f=4z, (16) 
\Xv Xv fiv) - Xj \v'J 



dx 
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where 

f{v) = v- v-^h{v) =2v- a{-f - l)v~^ - (a + 1), (17) 
with h as in (llip and a as in ([8]), or, equivalently, 



/(f-)=2^-h-l)(^J^^j (18) 

Eigenvalues of (jl2p correspond to nontrivial solutions W for which 
the boundary conditions W{±oo) = are satisfied. Because A{x, A) 
as a function of v is asymptotically constant in x, the behavior near 
X = iboo of solutions of (jl6p is governed by the limiting constant- 
coefficient systems 

W' = A±{\)W, A±{X) := A{±oo, A), (19) 

from which we readily find on the (nonstable) domain 5ReA > 0, 
A 7^ of interest that there is a one-dimensional unstable manifold 
(x) of solutions decaying at x = — oo and a two-dimensional stable 
manifold W2~{x) A W^{x) of solutions decaying at x = +oo, each of 
which may be chosen analytically in A. With additional care, these 
may be extended analytically to the whole set 5ReA > 0, i.e., to A = 
[13] . Defining the Evans function D associated with operator L as 
the analytic function 

D{X) := det{W^W+W+)^,=o, (20) 

we find that eigenvalues of L correspond to zeroes of D both in loca- 
tion and multiplicity; see, e.g., [H[T3 l [3 H H6] for further details. 

Equivalently, following |361[2] . we may express the Evans function 

as 

DW = {Wt-W^)^^^,, (21) 

where W^{x) spans the one-dimensional unstable manifold of solu- 
tions decaying at x = +oo (necessarily orthogonal to the span of 
W2{x) and Wr^{x)) of the adjoint eigenvalue ODE 

W' = -A{x,XyW. (22) 
The simpler representation (j2ip is the one that we shall use here. 



3. Description of the main results 



We can now state precisely our main results. 
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3.1. Limiting equations 

Under the strategic rescaling ([5]), both profile and eigenvalues equa- 
tions converge pointwise as — > to limiting equations at f+ = 0. 
The limiting profile equation (the limit of ([7])) is evidently 



v' = v{v - 1), (23) 



with explicit solution 



^ , , 1 - tanh(x/2) 

Mx) = 2 ' ^ ^ 

while the limiting eigenvalue system (the limit of (jl6p ) is 

W = A°{x,X)W, (25) 

/OA 1 \ 
^°(x,A) =00 1 , (26) 

\Xvo Xvo fo{vo) - Xj 



where 



Mvo) = 2vo-l = - tanh(x/2). (27) 



Indeed, this convergence is uniform on any interval vq > e > 0, or, 
equivalently, x < L, for L any positive constant, where the sequence is 
therefore a regular perturbation of its limit. We will call x G (— oo, L] 
the "regular region" or "regular side" . For — > on the other hand, 
or X — > oo, the limit is less well-behaved, as may be seen by the fact 
that df/dv ~ as {) — > f+, a consequence of the appearance of 
(^) in the expression (fTSi) for /. Similarly, in contrast to v, A{x,X) 
does not converge to A+(A) as x — > -|-oo with uniform exponential 
rate independent of f+, 7, but rather as Cv~^e~^/^ . We call x S 
[L, +00) therefore the "singular region " or "singular side" . (This is 
not a singular perturbation in the usual sense but a weaker type of 
singularity, at least as we have framed the problem here.) 



3.2. Limiting Evans function 

We should now like to define a limiting Evans function following the 
asymptotic Evans function framework introduced in [37J and establish 
convergence to this function in the u+ ^ limit, thus reducing the 
stability problem as f+ — > to the study of the (fixed) limiting 
Evans function. However, we face certain difficulties due to the (mild) 
singularity of the limit, as can be seen even at the first step of defining 
an Evans function for the limiting system. 
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For, the limiting coefficient matrix 



/OA 1 

(+00, A) =00 1 

Vo -1 - A 



) 



(28) 



is nonhyperbolic (in ODE sense) for all A, having eigenvalues 0, 0, — 1— 
A; in particular, the stable manifold drops to dimension one in the 
limit f+ 0. Thus, the subspace in which W2 and should be 
initialized at x = +00 is not self-determined by (f28]) . but must be 
deduced by a careful study of the double limit f+ — > 0, x — > +00. 
But, the computation 



lim ^(+00, A) =00 1 / Al{X) = lim lim A{x,X) 



shows that these limits do not commute, except in the special case 
7 = 1 already treated in [34J by other methods. 

The rigorous treatment of this issue is the main work of the paper. 
However, the end result can be easily motivated on heuristic grounds. 
A study of lim^i^^o ^(+00, A) on the set SReA > of interest reveals 
that eigenmodes decouple into a single "fast" (stable subspace) de- 
caying mode (*, *, 1)-^ associated with eigenvalue —7 — A of strictly 
negative real part and a two-dimensional (center) subspace of neutral 



only a single genuine eigenvector (1,0,0)"^. For small, therefore, 
j4_|_(A) has also a single fast, decaying, eigenmode with eigenvalue 
near —7 — 1 and two slow eigenmodes with eigenvalues near zero, 
one decaying and the other growing (recall. Section 12.41 that the sta- 
ble subspace of has dimension two for "SieX > 0, A 7^ and the 
unstable subspace dimension one). 

Focusing on the single slow decaying eigenvector of A^, and con- 
sidering its limiting behavior as f+ 0, we see immediately that it 
must converge in direction to ±(1,0,0)^. For, the sequence of direc- 
tion vectors, since continuously varying and restricted to a compact 
set, has a nonempty, connected set of accumulation points, and these 
must be eigenvectors of lim^^^o ^+ with eigenvalues near zero. Since 
±(1,0,0)"^ are the unique candidates, we obtain the result. Indeed, 
both growing and decaying slow eigenmodes must converge to this 
common direction, making the limiting analysis trivial. 

The same argument shows that ±(1, 0, O)"'" is the limiting direction 
of the slow stable eigenmode of A^{x, X) as x — > -|-oo, that is, in the 
alternate limit lima;_>oo lini^+^o ^(2;, A). That is, := (1,0,0)-^ is 
the common limit of the slow decaying eigenmode in either of the 




(29) 



modes (r, 0) 



associated with Jordan block 
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two alternative limits lim^_|^-^o ^+ and lim»;4.-^o it thus seems a 
reasonable choice to use this limiting slow direction to define an Evans 
function for the limiting system (j26l) . On the other hand, the stable 
eigenmode 

^3 := {a-\\/a + l),a-\Vf, 

a = —1 — A, of A^j^ is forced on us by the system itself, independent 
of the limiting process. 

Combining these two observations, we require that solutions 
and of the limiting eigenvalue system (|26p lie asympotically in 
directions V2 and V3, respectively, thus determining a limiting, or 
"reduced" Evans function 

D\\) := det«-H^°+W^3°+)|,=o, (30) 

or alternatively 

^°(A) = (^^^{'+.<-)|^^„, (31) 

with W^~^ defined analogously as a solution of the adjoint limiting 
system lying asymptotically at x = +00 in direction 

Vl:=(0,l,a-i)^ = (0,l,(l + A)-i)^ (32) 

orthogonal to the span of V2 and V3, where "~" denotes complex 
conjugate. (The prescription of in the regular region is straight- 
forward: it must lie on the one-dimensional unstable manifold of 
as in the > case.) 



3.3. Physical interpretation 

Alternatively, the limiting equations may be derived by taking a for- 
mal limit as w+ ^ of the rescaled equations ([6]), recalling that 
a ~ w^, to obtain a limiting evolution equation 

Vt + Vx - = 0, 

(Ux\ (33) 

corresponding to a pressure-less gas, or 7 = 0, then deriving profile 
and eigenvalue equations from (j33p in the usual way. This gives some 
additional insight on behavior, of which we make important mathe- 
matical use in Appendix[Bj Physically, it has the interpretation that, 
in the high-Mach number limit — > 0, effects of pressure are con- 
centrated near x = +00 on the infinite-density side, as encoded in 
the special upstream boundary condition {u,u' ,v,v') — > c(l, 0,0,0) 
as X — > -|-oo for the integrated eigenvalue equation, which may be seen 
to be equivalent to the conditions imposed on W^' in the previous 
subsection. 
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3.4- Analytical results 

Defining as in (|30p -(|3i p . we have the following main theorems. 

Theorem 1. For A in any compact subset of ?R.e\ > 0, -D(A) con- 
verges uniformly to D^(X) as Vj^ ^0. 

Corollary 1. For any compact subset A of 5ReA > 0, D is nonvan- 
ishing on A for sufficiently small if is nonvanishing on A, and 
is nonvanishing on the interior of A only if is nonvanishing there. 

Proof. Standard properties of uniform limit of analytic functions. 

Corollary 2. Isentropic Navier-Stokes shocks are stable in the high- 
Mach number limit f + ^ i/ is nonvanishing on the wedge 



and only if D is nonvanishing on the interior of A. 

Proof. Corollary [1] together with Proposition [2j 

Remark 1. Likewise, on any compact subset of ?R.eX > 0, \D\ is uni- 
formly bounded from zero for sufficiently small {M sufficiently 
large) if and only if \D^\ is uniformly bounded from zero. Thus, isen- 
tropic Navier-Stokes shocks are "uniformly stable" for sufficiently 
small v^, in the sense that |D| is bounded from below independent 
of V-^-, if and only if is nonvanishing on A as defined in ()34p . 

The following result completes our abstract stability analysis. The 
proof, given in Appendix [Bl is by an energy estimate analogous to 
that of [MlEI 

Proposition 3. The limiting Evans function is nonzero on 5ReA > 



Corollary 3. For any 7 > 1, isentropic Navier-Stokes shocks are 
stable for Mach number M sufficiently large (equivalently, suffi- 
ciently small). 



^ Stability for 7 = 1, proved in [34) . already implies nonvanishing of D'^ outside 
the imaginary interval [—i^3/2, +iy^3/2], by Corollary [2] 





(34) 



0. 
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3.5. Numerical computations 

Unfortunately, the energy estimate used to establish [3l though math- 
ematically elegant, yields only the stated, abstract result and not 
quantitative estimates. A simpler and more general approach, that 
does yield quantitative information, is to compute the reduced Evans 
function numerically. We carry out this by-now routine numerical 
computation using the methods of [2]. Specifically, we map a semi- 
circle d{{^eX > 0} n {|A| < 10}) enclosing yl for 7 G [1,3] by 
and compute the winding number of its image about the origin to 
determine the number of zeroes of within the semicircle, and thus 
within A. For details of the numerical algorithm, see [21 111^121] . 

The result is displayed in Figure [U clearly indicating stability. 
More precisely, the minimum of |-D| on the semicircle is found to be ~ 
0.2433. Together with Rouches Theorem, this gives explicit bounds on 
the size of the Mach number for which shocks are stable, as displayed 
in Tabled Section E specifically, M > 50 for 7 G [1,2], M > 2,400 
for 7 G [2, 2.5], and M > 13, 000 for 7 G [2.5, 3], all corresponding ap- 
proximately to u+ = 10"'^. In Figure [21 we superimpose on the image 
of the semicircle by its (again numerically computed) image by 
the full Evans function Z), for a monatomic gas 7 = 1.66... at succes- 
sively higher Mach numbers u+ =le-l , le-2 , le-3, le-4, le-5 , le-6, 
graphically demonstrating the convergence of D to as V-\- ap- 
proaches zero. 

Indeed, we can see a great deal more from Figure [2j For, note that 
the displayed contours are, to the scale visible by eye, "monotone" in 
1)+, or nested, one within the other (they do not appear so at smaller 
scales). Thus, lower-Mach number contours are essentially "trapped" 
within higher-Mach number contours, with the worst-case, outmost 
contour corresponding to the limiting Evans function D^. From this 
observation, we may conclude with confidence stability down to the 
smallest value M ~ 5.5 displayed in the figure, corresponding (see 
Table [2]) to f+ = 10~^. That is, a great deal of topological informa- 
tion is encoded in the analytic family of Evans functions indexed by 
from which stability may be deduced almost by inspection. Be- 
havior for other 7 G [0, 3] is similar. See, for example, the case 7 = 3 
displayed in Figure |3l which is virtually identical to Figure [2FI 

Such topological information does not seem to be available from 
other methods of investigating stability such as direct discretation 
of the linearized operator about the wave [26] or studies based on 
linearized time-evolution or power methods This represents in 

our view a significant advantage of the Evans function formulation. 



In particular, Figure [3] indicates stability down to 11+ = 10 ^, or Mach number 
~ 20, from which we may conclude unconditional stability on the whole range 

7e [1,3] of [g. 
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Fig. 1. The image of the semi-circular contour via the Evans function for the 
reduced system. Note that the winding number of this graph is zero. Hence, there 
are no unstable eigenvalues in the semi-circle. 



Remark 2. Recall that the Evans function is not determined uniquely, 
but only up to a nonvanishing analytic factor [I],ll3j . The simple con- 
tour structure in Figure [2] is thus partly due to a favorable choice 
of D induced by the initialization at ±00 by Kato's ODE [24], as 
described in |li y221 [2]. A canonical algorithm for tracking bases of 
evolving subspaces, this in some sense minimizes "action"; see |20j 
for further discussion. 

Remark 3. Note that the limiting equations, and the limiting Evans 
function are both independent of 7. To study high-Mach number 
stability for 7 > 3, therefore, requires only to examine on suc- 
cessively larger semicircles. Thus, our methods in combination with 
the those of [2], allow us to determine stability in principle over any 
bounded interval in 7, for 7 > 1 and for all Mach numbers M > 1. 

Remark 4- As Figure [2] suggests, an alternative method for determin- 
ing stability, without reference to D^, is to compute the full Evans 
function for sufficiently high Mach number. That is, nonvanishing of 
D^, and thus stability of sufficiently high Mach number shocks for 
7 G [1,3], can already be concluded by large-but-finite Mach number 
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Fig. 2. Convergence to the limiting Evans function as — > for a monatomic 
gas, 7 = 1.66.... The contours depicted, going from inner to outer, are images of 
the semicircle under D for v+ =le-l , le-2, le-3, le-4, le-5, le-6. The outermost 
contour is the image under 73°, which is nearly indistinguishable from the image 
for v+ =le-6. 

study of [2] together with the fact that a hmit D — > exists (see 
also Remark d]). 



3.6. Conclusions 

The analytical result of Corollary [3] guarantees stability for 7 > 1, 
M sufficiently large. For 7 G [1,2.5], our numerical results indicate 
stability for M > 2, 500 by a crude Rouche bound, and indeed much 
lower if further structure is taken into account. Together with the 
small and intermediate Mach number studies of [Mill] for M < 3, 000, 
this yields unconditional stability of isentropic Navier-Stokes shocks 
for 7 G [1,2.5] and M > 1. There is no inherent restriction to 7 G 
[1,2.5]; as discussed in Remark [3l numerical computations can be 
carried out for any value of 7 to determine stability (or instability) 
for all M > 1. Indeed, our method of analysis indicates that the 
large-7 limit is quite analogous to the high-Mach number limit f + — > 
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Fig. 3. Convergence to the limiting Evans function as u+ — > for 7 = 1. The 
contours depicted, going from inner to outer, are images of the semicircle under D 
for v+ =le-l , le-2, le-3 , le-4, le-5 , le-6. The outermost contour is the image 
under _D°, which is nearly indistinguishable from the image for v+ =le-6. 



0, suggesting the possibility to establish still more general results 
encompassing all 7 > 1, M > 1. 

Our numerical results reveal also an unexpected "universality" of 
behavior in the high-Mach number regime, beyond just convergence 
to the limiting system. Namely, we see (cf. Figures [2] and H]) that 
behavior for a given f+ is virtually independent of the value of 7. 
This also indicates that u + and not M is the more useful measure of 
shock strength in this regime. 



4. Boundary-layer analysis 

We now carry out the main work of the paper, analyzing the flow of 
(|16|) in the singular region. Our starting point is the observation that 

/OA 1 \ 

A{x,\) = 1 (35) 
\\v Xv f{v) - X) 
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Fig. 4. Convergence to the limiting Evans function as u+ — > for 7 = 3. The 
contours depicted, going from inner to outer, are images of the semicircle under D 
for vj^ =le-l , le-2, le-3 , le-4, le-5 , le-6. The outermost contour is the image 
under _D°, which is nearly indistinguishable from the image for v+ =le-6. 

is approximately block upper-triangular for v sufficiently small, with 

diagonal blocks and {f{v) — A) that are uniformly spectrally 

separated on 5?eA > 0, as follows by 

f{v) <2v -l<-l/2. (36) 

We exploit this structure by a judicious coordinate change convert- 
ing (jl6p to a system in exact upper triangular form, for which the 
decoupled "slow" upper lefthand 2x2 block undergoes a regular per- 
turbation that can be analyzed by standard tools introduced in |37| . 
Meanwhile, the fast, lower righthand 1x1 block, since scalar, may 
be solved exactly. 

The global structure of this argument loosely follows the general 
strategy of [37j of first decoupling fast and slow modes, then treat- 
ing slow modes by regular perturbation methods. However, there are 
interesting departures that may be of use in other degenerate situa- 
tions. First, we only partially decouple the system, to block-triangular 
rather than block-diagonal form as in more standard cases, and sec- 
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ond, we introduce a more stable method of block-reduction taking ac- 
count of usually negligible derivative terms in the definition of block- 
triangularizing transformations, which, if ignored, would in this case 
lead to unacceptably large errors. 



4.1. Preliminary transformation 

We first block upper-triangularize by a static (constant) coordinate 
transformation the limiting matrix 

/OA 1 \ 

A+ = ^(+00, A) = 1 (37) 
\Xv+ Xv+ fiv+)- Xj 

at X = -|-cxD using special block lower-triangular transformations 

^+ ■= {xvi9+ 1) ' ^+ •= = {-xv+e+ 1) ' 

where / denotes the 2x2 identity matrix and 0+ G C"*^^^ is a 1 x 2 
row vector. 

Lemma 1. On any compact subset of ?R.eX > 0, for each v+ > 
sufficiently small, there exists a unique 9^ = ^+(f+,A) such that 
:= L+^4_i2-|- is upper block-triangular, 

where J = ^''^^ ~ ' satisfying a uniform bound 

\9+\<C. (40) 

Proof. Setting the 2 — 1 block of to zero, we obtain the matrix 
equation 

9+{aI - AJ) = -1^ + XV+9+19+, 
where a = f{v+) — A, or, equivalently, the fixed-point equation 

9+ = {aI-Xjy^(-l^ + XV+9+19+). (41) 



By det(al — XJ) = / 0, {al — XJ)^^ is uniformly bounded on 
compact subsets of 5ReA > (indeed, it is uniformly bounded on 
all of 3?eA > 0), whence, for |A| bounded and sufficiently small, 
there exists a unique solution by the Contraction Mapping Theorem, 
which, moreover, satisfies (1401). 
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4-2. Dynamic triangularization 
Defining now Y := L^W and 
A{x, A) = L+A{x, X)R+{x, A) 

X{v - v+)l^ - Xv+if{v) - f{v+))e+ fiv) - A - 

(42) 

we have converted (I16p to an asymptotically block upper-triangular 
system 

Y' = i(x, A)y, (43) 

with j4+ = A{+oo, A) as in ()39p . Our next step is to choose a dynamic 
transformation of the same form 



converting (j43p to an exactly block upper-triangular system, with 
uniformly exponentially decaying at x = -|-oo: that is, a regular 
perturbation of the identity. 

Lemma 2. On any compact subset of SReA > 0, for L sufficiently 
large and each > sufficiently small, there exists a unique = 
0+{x, A, W-|_) such that A := LA{x, X)R+L' R is upper block-triangular, 

^"V f{v)-x-xe+\-0\)' 

and 0{L) = 0, satisfying a uniform bound 

|0(x,A,u+)| < Ce"''^, r]>0,x>L, (46) 
independent of the choice of L, v+. 

Proof. Setting the 2 — 1 block of A to zero and computing 

0\ // 0\ / 0' 



' -0' Oj \0 Ij \-0' 0, 
we obtain the matrix equation 

0' -0{aI-X{J + v+19+))=C + 01il0, (47) 
where the forcing term 

^ := -X{v - v+)l^ + Xv+{fiv) - fiv+))e+ 
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by derivative estimate df /dv < Cv ^ together with the Mean Value 
Theorem is uniformly exponentially decaying: 

Id < C|t)-w+| < Cae"''^, ?7>0. (48) 

Initializing 0{L) = 0, we obtain by Duhamel's Principle/ Variation 
of Constants the representation (supressing the argument A) 

e{x) = 1^ Sy^^iC + 01O){y) dy, (49) 

where S^^^ is the solution operator for the homogeneous equation 
O' -0{al- A(J + V+19+)) = 0, 

or, explicitly, 

gy^x ^ ^fy a{y)dy ^-X(J+v+l_e+)ix-y) ^ 



For |A| bounded and u+ sufficiently small, we have by matrix per- 
turbation theory that the eigenvalues of — A(J + f+10+) are small 
and the entries are bounded, hence 

for z > 0. Recalling the uniform spectral gap 5Rea = f{v) — 5ReA < 
— 1/2 for KeA > 0, we thus have 

ISy^'^l < C7e^(^-^) (50) 
for some C, > 0. Combining ([18|) and ([50]) . we obtain 



y 



C7 e-(r?/2)x_ 



(51) 



Defining 0{x) =: 6(x)e and recalling (j49p we thus have 

0{x) = f + e^''/^)^ sy^^'e-'^yeieiy) dy, (52) 

where / := e^^^"^^^ Sy^^C{y) dy is uniformly bounded, |/| < C3, 
and e'-''/^-*^ Sy~*^e~^y91\.9{y) dy is contractive with arbitrarily small 
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contraction constant e > in L°°[L, +00) for 16*1 < 2C3 for L suffi- 
ciently large, by the calculation 

Jl Jl 

^{ri/2)x I' i^^-ri(x-y)^-r,y 
(j^-{^/2){x~y) 



< 



L 



1 — "2II00 max Ipjlloo 
j 



< e 



-{ri/2)L 



L 



2 00 max (7 J 00 
j 



= Cge-^''/^)^!!^! - ^2||oomax||6*j|U. 

j 

It follows by the Contraction Mapping Principle that there exists a 
unique solution 9 of fixed point equation (j52p with \6{x)\ < 2C3 for 
X > L, or, equivalently (redefining the unspecified constant r/), (|16]) . 

Remark 5. The above calculation is the most delicate part of the anal- 
ysis, and the main technical point of the paper. The interested reader 
may verify that a "quasi-static" transformation treating term O' in 
(I47p as an error, as is typically used in situations of slowly-varying co- 
efficients (see for example [3HI37]). would lead to unnaceptable errors 
of magnitude 

Oi\f{v)y\\v - v+\) = 0{\df/dv)\\v - v+\) = 0{\v\-'\v - v+\). 

One may think of the exact ODE solution ()44p as "averaging" the 
effects of rap idly- varying coefficients by integration of (j47p . 



4.3. Fast/Slow dynamics 

Making now the further change of coordinates 

Z = LY 

and computing 

(LY)' = LY' + L'Y = {LA+ + L')Y, 

= {LA+R + L'R)Z, 

we find that we have converted (j43p to a block-triangular system 

f{v) - A - \v+e+\ - 01, 

(53) 
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related to the original eigenvalue system (I16|) by 

W = LZ, R := R+R = , L := R-' = {j^ , (54) 

where 

e = e + \v+e+. (55) 

Since it is triangular, (j53p may be solved completely if we can solve 
the component systems associated with its diagonal blocks. The fast 
system 

z' = (f{v) - A - \v+e+\ - e\)z 



associated to the lower righthand block features rapidly-varying co- 
efficients. However, because it is scalar, it can be solved explicitly by 
exponentiation. 
The slow system 

z' = (\{J + v+\e+) + \e]z (56) 



associated to the upper lefthand block, on the other hand, by (l46l) . 
is an exponentially decaying perturbation of a constant-coefficient 
system 

z' = \{J + v+\e+)z (57) 

that can be explicitly solved by exponentiation, and thus can be well- 
estimated by comparison with ()57p . A rigorous version of this state- 
ment is given by the conjugation lemma of [35j: 

Proposition 4 ([35]). Let M{x,X) = M+(A) + 6'(x,A), with M+ 
continuous in A and |6'(a;,A)| < Ce~^^ , for A in some compact 
set A. Then, there exists a globally invertible matrix P{x,X) = I + 
Q{x,X) such that the coordinate change z = Pv converts the variable- 
coefficient ODE z' = M{x,X)z to a constant- coefficient equation 

v' = M+(X)v, 

satisfying for any L, < fj < rj a uniform bound 

|Q(x,A)| < C(L,r/,r/,max|(M+)ij|,dimAf+)e-'^^ for x > L. (58) 

Proof. See [35y46j. or Appendix ICl 

By Proposition m the solution operator for (j56p is given by 

P{y, A)e^(^+^+i'^+(^'^+))(^-s^)p(x, A)^\ (59) 

where P is a uniformly small perturbation of the identity for x > L 
and L > sufficiently large. 
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5. Proof of the main theorem 

With these preparations, we turn now to the proof of the main the- 
orem. 



5.1. Boundary estimate 

We begin by establishing the fohowing key estimates on W^{L), that 

is, the value of the dual mode appearing in ([2T]) at the boundary 
X = L between regular and singular regions. 

Lemma 3. For A on any compact subset of 3fteA > 0, and L > 
sufficiently large, with normalized as in \1S\ 31 W^. 

-^il < Ce"^^ (60) 
as u+ 0, uniformly in X, where C, r] > are independent of L and 

Vi := (0,l,(l + A)-i)^ 
is the limiting direction vector (I32p appearing in the definition of . 
Corollary 4. Under the hypotheses of Lemma\^ 

\W^+{L)-Vi\<Ce-'>^ (61) 

and 

\W+{L) - W^+{L)\ < Ce-'i^ (62) 

as 0, uniformly in A, where C, r] > are independent of L 

and is the solution of the limiting adjoint eigenvalue system 

appearing m definition of . 

Proof (Proof of Lemma\B^. Making the coordinate-change 

Z := R*W, (63) 

R as in (f54l) . reduces the adjoint equation W' = —A*W to block 
lower-triangular form, 

Z' = -A*Z 

-xj^ + {Xv+i9+ - ley _ _ \ 

-1^ -f{v) + X + Xv+{9+l)*-{0irJ ' 

(64) 

with " " denoting complex conjugate. 
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Denoting by a suitably normalized element of the one-dimensional 
(slow) stable subspace of —A*, we find, similarly as in the discussion 
of Section [3T2] that, without loss of generality, 

y+-(0,l,(7 + A)-i)^ (65) 

as t;+ ^0, while the associated eigenvalue flf —>■ 0, uniformly for A 
on an compact subset of 5ReA > 0. The dual mode = R*W^ is 
uniquely determined by the property that it is asymptotic as x — > 

+00 to the corresponding constant-coefficient solution e'^i (the 
standard normalization of p^[371 [2]). 

By lower block-triangular form (j64[) . the equations for the slow 
variable z'^ := (Zi, Z2) decouples as a slow system 

z' = -(^X{J + v+19+) + ieyz (66) 

dual to (156]) . with solution operator 

P*{x, A)-ie-^('^+^+^^+)*)(^-2^)p(y, A)* (67) 
dual to (f59]) . i.e. (fixing y = L, say), solutions of general form 

z(A, x) = P*{x, A)-ie-^(-^+''+^^+)*)(^-^)t), (68) 

■u G arbitrary. 
Denoting by 

Z+(L) := R*W^{L), 

therefore, the unique (up to constant factor) decaying solution at 
-|-oo, and vf := {{V{^)i, {V^)2)'^ , we thus have evidently 

z+{x, A) = P*{x, A)-ie-^(-^+''+^^+)*)^t}+, 

which, as 0, is uniformly bounded by 

|5+(x,A)| < Ce^^ (69) 

for arbitrarily small e > and, by (j65p . converges for x less than or 
equal to any fixed X simply to 

\im zt{x,\) = P*{x,X)-\0,lf. (70) 
1;+— >o 

Defining by q := {Z^)^ the fast coordinate of Z^ , we have, by 
q' + {f{v) - A - (Xv+e+l + 01)*) g = 
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whence, by Duhamel's principle, any decaying solution is given by 

J X 

where 

a{y, A, v+) := - (/(?)) - A - {Xv+d+1 + 01)*) . 

RecalUng, for ReX > 0, that ^ea > 1/2, combining (IM]) and UiTUtf . 
and noting that a converges uniformly on y < Y as ^ for any 
y > to 

ao{y,X) :=-/o(t)) + A + (0oir 

= {i + x) + o{e-''y) 

we obtain by the Lebesgue Dominated Convergence Theorem that 
q{L, A) ^ y e^y M-'^)dz^T^Q^ ^^T 



L 

= (l + A)-i(l + 0(e-''^)). 

Recalling, finally, (j70p . and the fact that 

|P-/d|(L,A), |i?-/d|(L,A) < Ce-'?^ 

for f+ sufficiently small, we obtain (j60p as claimed. 

Proof (Proof of Corollary^. Applying Proposition U] to the limiting 
adjoint system 

W' = -{A^yW = I -A _\W + 0{e-'>^)W, 




we find that, up to an Id + 0{e coordinate change, 
is given by the exact solution W = Vi of the limiting, constant- 
coefficient system 



w' = -{A'^yw = 

This yields immediately ()6ip . which, together with (I60p . yields ([62 
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Remark 6. Noting that ()6ip is sharp, we see from (I62|) that the error 
between W^{L) and W^'^{L) is already within the error tolerance 
of the numerical scheme used to approximate , in which W^^ is 
initialized at x = L with approximate value Vi [TTllSTl ll]. Thus, so 
long as the flow on the regular region x < L well-approximates the 
exact limiting flow as f+ 0, we can expect convergence of D to 
based on the known convergence of the numerical approximation 
scheme. 



5.2. Convergence to 

As hinted by Remark [H the rest of our analysis is standard if not 
entirely routine. 

Lemma 4. On x < L for any fixed L > 0, there exists a coordinate- 
change W = TZ conjugating (jl6p to the limiting equations (|26p . 
T = r(x. A, satisfying a uniform bound 

\T - Id\ < C{L)v+ (71) 

for all u-i- > sufficiently small. 

Proof. For x G (— oo,0], this is a consequence of the Convergence 
Lemma of [37j, a variation on Proposition HI together with uniform 
convergence of the profile and eigenvalue equations. For x G [0, L], it 
is essentially continuous dependence; more precisely, observing that 
\A - < Ci{L)v+ for x G [0,L], setting S := T - Id, and writing 
the homological equation expressing conjugacy of (jl6p and (j26p . we 
obtain 

S'-iAS-SA'^) = (A-A^), 

which, considered as an inhomogeneous linear matrix- valued equa- 
tion, yields an exponential growth bound 

S{x) < e^^(5(0) + C-'^Ci{L)v+) 
for some C > 0, giving the result. 

Proof (Proof of Theorem\^. Lemma lU together with convergence as 
f + ^ of the unstable subspace of A_ to the unstable subspace of 
A^_ at the same rate 0(^^+) (as follows by spectral separation of the 
unstable eigenvalue of AP and standard matrix perturbation theory) 
yields 

|VFf (0, A) - VF°-(0,A)| < C{L)v+. (72) 
Likewise, Lemma H] gives 

|H^+(0,A) -M^o+(0,A)| < C{L)v+\Wt{QA)\ 

+ |5o^^°||Ty+(L,A)-H^r(L,A)|, 
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where Sq~^^ denotes the solution operator of the hmiting adjoint 
eigenvalue equation W' = —{A^)*W. Applying Proposition H] to the 
limiting system, we obtain 

1^0^° I < Cae-^'i^ < C2L\X\ 

by direct computation of e~^+^, where C2 is independent of L > 0. 
Together with ([62]) and ([73]) . this gives 

|VF+(0,A) -T^°+(0,A)| < CiL)v+\W+iO,X)\+ 11X10206-"^, 
hence, for |A| bounded, 

|VF+(0,A) - W^+iO,X)\ < Cs{L)v^\W^+{0,X)\+LC,e-'^'^ 

< C5iL)v+ + LC^e-^^^ . 

Taking first L — > co and then — > 0, we obtain therefore conver- 
gence of W^{0, A) and W^{0, A) to W^+{0, A) and W^+{0, A), yield- 
ing the result by definitions ([21|) and ([3ip . 



6. Numerical convergence 

Having established analytically convergence of D to D° as M — > 00, 
we turn finally to numerics to obtain quantitative information yield- 
ing a concrete stability threshold. Specifically, for fixed 7, we com- 
pute the "Rouche bound" v+ at which the maximum relative error 
\D — D^\/\D^\ over the semicircular contour d{^eX > 0, |A| < 10} 
around which we perform our winding number calculations becomes 
1/2. Recall that Rouche's Theorem guarantees for relative error < 1 
that the winding number of D is equal to the winding number of D^, 
which we have shown to be zero, hence we may conservatively con- 
clude stability for f+ less than or equal to this bound, or M greater 
than or equal to the corresponding Mach number. Computations are 
performed using the algorithm of [2]; results are displayed in Table 

m 

More detailed results are displayed for the monatomic gas case, 
7 = 1.66..., in TableO Results are similar for other 7 E [1, 3], as may 
be seen by comparing Figures [21 [3l and [4l 

Prom the quantitative gap and conjugation estimates given in Ap- 
pendix [C] from which follows also a quantitative version of the Con- 
vergence Lemma of [37], one could in principle establish quantita- 
tive convergence rates for \D — D^\, by tracking constants carefully 
through the estimates of the previous sections. Indeed, one could do 
much better than the rather crude bounds stated for the general case 
by taking into account the eigenstructure of the actual matrices A± , 
appearing in our analysis. That is, there are contained in our 
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Relative 


Mach 


7 




Error 


Number 


3.0 


1.27e-3 


.5009 


12765 


2.5 


1.36e-3 


.5006 


2423 


2.0 


1.49e-3 


.5001 


474 


1.5 


1.75e-3 


.4999 


95.5 


1.0 


2.8e-3 


.4995 


18.9 



Table 1. Rouche bounds for various 7. 





Mach 


Relative 


Absolute 


v+ 


Number 


Difference 


Difference 


1.0(-6) 


7.71(4) 


0.1221 


0.0601 


1.0(-5) 


1.13(4) 


0.1236 


0.1445 


1.0(-4) 


1.64(3) 


0.1487 


0.4714 


1.0(-3) 


2.44(2) 


0.4098 


1.3464 


1.0(-2) 


36.1 


0.9046 


2.8253 


l.O(-l) 


5.50 


1.2386 


3.8688 



Table 2. Maximum relative and absolute differences between D and D , for 
7 = 1.66... and A on the semicircle of radius 10. 



analysis, as in the study of [2j, all of the ingredients needed for a nu- 
merical proof. Given the fundamental nature of the problem studied, 
this would be a very interesting program to carry out. 



7. Discussion and open problems 

Besides long-time stability, our results have application also to ex- 
istence of shock layers in the small- viscosity limit, which likewise 
reduces to the question of stability of the Evans function [38..17j . 
Indeed, spectral stability has been a key missing piece in several di- 
rections [45','46]. Our methods should have application also to spectral 
stability of large-amplitude noncharacteristic boundary layers, com- 
pleting the investigations of [iTlflGlISS] . It may be hoped that they 
will extend also to full gas dynamics and multi-dimensions, two im- 
portant directions for further investigation. As discussed in the text, 
the problems of numerical proof and of stability in the large-7 limit 
are two other interesting directions for further study. 

More speculatively, our results suggest the possibility of a large- 
variation version of the results obtained by quite different methods 
in [7] on general viscous solutions (including not only noninteracting 
shocks, but shocks, rarefactions, and their interactions), and, through 
the physical insight provided into the high-Mach number limit, per- 
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haps even a hint toward possible methods of analysis. This would be 
an extremely interesting direction for further investigation. 



A. Proofs of Preliminary Estimates 

Proof (Proof of Proposition [7]j. Existence and monotonicity follow 
trivially by the fact that ([7]) is a scalar first-order ODE with convex 
righthand side. Exponential convergence as a; — > +00 follows, for 
example, by the computation 

)-(HT)(e)^-0) 



yielding 

H{v,v+) 
f — 7 < < V — (I — v+) 

V — v+ 

by the elementary estimate 1 < < 7 for < x < 1 . Convergence 
as X — > —00 follows by a similar, but simpler computation; see p]. 

Lemma 5. The following identity holds for 5ReA > 0.' 

(5ie(A) + |3m(A)|) / M'^ -\ I + [ 

Jm. ^ Jm. Jr 

<^ [ [ y\u'\\u\. (75) 

Jr 

Proof. We multiply (jl2b[) by vu and integrate along x. This yields 

\ /" 'I |2 I /" - , /" I /|2 f Hv) 
A / v\u\ + VUU+ \u\ = — — V u. 

Jr Jr Jr Jr v'^ 

We get ([75]) by taking the real and imaginary parts and adding them 
together, and noting that |3^e(z)| + |9m(z)| < \/2|z|. □ 

Lemma 6. The following identity holds for JReA > 0.' 

/ = 25Re(A)2 / |t;|2 + 5Re(A) / ^ + i / 
Jr Jr Jr v 2 

(76) 



hjv) 



+ 



07 



,,'|2 
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Proof. We multiply (|12b|) by v' and integrate along x. This yields 



A / uv' + / uv 



Kv) I /|2 



1 

— U V 
V 



1 / \ / I II \ -I 
-{AV + V jV . 

R V 



Using (|12ap on the right-hand side, integrating by parts, and taking 
the real part gives 



A / uv' + / u'v 



h{v) _^ 



{)7+l 2£'2 

The right hand side can be rewritten as 



+ 3f?e(A) [ 

Jr V 



A / Mu' + / u'v 



h{v) ^ a-f 



{;7+l 1)7+1 

Now we manipulate the left-hand side. Note that 



v'\^ + ^eiX) [ 

Jr V 
(77) 



A / uv + I u'v' = (A -I- A) / uv' — / u{\v' + v") 

Jr Jr 

= -25Re(A) [ u'v - [ uu" 
Jr Jr 

= -23f?e(A) / iXv + v')v+ [ |n'p. 



Hence, by taking the real part we get 



5Re 



A / uv' + / u'v 



-25Re(A)" / \v 



This combines with (|77p to give (j76p . □ 
Lemma 7. For /i('u) as in (jlip . «;e /laue 



sup 

i) 



7- 



<7- 



(78) 



Proof. Defining 

giv) := h{v)v-^ = -v + a{j - l)v~'^ + {a + 1), (79) 

we have g'{v) = — 1 — 07(7 — l)v~^~^ < for < u+ < t) < f_ = 1, 
hence the maximum oi g on v E [f +,?;_] is achieved at v = Vj^.. 
Substituting ([8]) into ([79]) and simplifying yields ([78|) . 
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Proof (Proof of Proposition Wjj. Using Young's inequality twice on 
right-hand side of (f75]) together with ([78l) . we get 



(3fie(A) + |3=m(A)|) / v\uY - - / Vxlnr + / \v: 



l\2 



< V2 



v'y 



V \\u\ + / v\u \\u 



< 



< 



^!7^ 



/|2 



v\u\ +e / v\u \ + — / v\u 



4e 



h{v) ^, 



V \ + e / \u 



'|2 



20 4e 



v\u\ 



Assuming that < e < 1 and 9 = {1 — e)/2, this simphfies to 



(5Re(A) + |9m(A)|) / v\u\^ + {I - e) / |n 



./|2 



< 



1 - e /■ /i(u) , , 



r + 



7^1 
29 4e 



Applying ([76]) yields 

(Ke(A) + |9m(A)|) / v\u\'^ < 
Jr 

or equivalently, 

(3fie(A) + |9m(A)|) < 



7 1 



1 - e 4e 



(47 - l)e - 1 



V u 



4e(l - e) 



Setting e = 1/(2^ + 1) gives ([M]). □ 



B. Nonvanishing of D 

As pointed out in Section |3.3] the limiting eigenvalue system ()25p . 
()26p . together with the limiting boundary conditions derived in Sec- 
tion 13.21 may be expressed equivalently as the integrated eigenvalue 
problem 

Xv + v' -u' = 0, (80a) 
Xu + u' -^^v' = (80b) 

V V 

corresponding to a pressureless gas, 7 = 0, with special boundary 
conditions 



(u, u', V, v'){-oo) = (0, 0, 0, 0), {u, u',v, v'){+oo) = (c, 0, 0, 0). (81) 
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Motivated by this observation, we establish stabihty of the hmit- 
ing system by a Matsumura-Nishihara-type spectral energy estimate 
exactly analogous to that used to prove stability for 7 = 1 in 



Proof (Proof of Proposition\3^. Multiplying (iSObl) by vu/{l — v) and 
integrating on some subinterval [a.h] C M, we obtain 



A 



1 -v' 



\u\ dx 



1 



-u'udx 



v'udx 



^ u"u 



l-v 



-dx. 



Integrating the third and fourth terms by parts yields 



l-v 



\u\ dx + 



- + 



+ 



vu + 



1 - v' 
u'u 



1 



l-v 

■b 



u'udx 



dx + v{Xv + v')dx 



l-v 



Taking the real part, we have 



rb 1^/12 



3f?e(A) J + (ix + J g{v)\ufdx + J j-^dx 

\2- 



vu + 



uu 1 



l-v 2 



V I 1 

+ 



where 



9{v) 



1 — V \1 — V 



+ 



\V\ 



(82) 



l-v 



l-v 



Note that 
d f 1 



{l-v)' 



v{v — 1) 



dx\l-v) {l-v)2 (l-v)2 {l-v)2 



l-v 



Thus, g{v) = and the third term on the right-hand side vanishes, 
leaving 



3f?e(A) 



u\ +\v\ ] dx + I z rux 



l-v 



l-v 



?fie{vu) + 



'Sie{u'u) \v 



21 



l-v 



We show next that the right-hand side goes to zero in the limit 
as a ^ —00 and b 00. By Proposition HJ the behavior of u, v near 
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iboo is governed by the limiting constant-coefficient systems W' = 
A'].{X)W, where W = {u,v,v')'^ and = A^{±oo, A). In particular, 
solutions W asymptotic to (1,0,0) at x = +00 decay exponentially 
in [u' ,v,v') and are bounded in coordinate u as x — > +00. Observing 
that 1 — V lasx— > +00, we thus see immediately that the 
boundary contribution at b vanishes as 6 — > +00. 

The situation at —00 is more delicate, since the denominator 1 — v 
of the second term goes to zero at rate as x —00, the rate 
of convergence of the limiting profile v. By inspection, the limiting 
coefficient matrix 

/OA 1 \ 
A^_= {0 1 , (83) 
\A A 1 - A/ 



has eigenvalues 



12 = -A, 



1 ± VI + 4A 



2 

hence for SReA > the unique decaying mode at a: = +00 is the 
unstable eigenvector corresponding to ^ = l+v^+M^ with growth 

Se(l±4±i?) = l^s,^TT4A>l/2. 

Thus, \u'\, \v'\, \v\ < Ce^^'^^^^^'^ as x ^ —00, e > 0, and in partic- 
ular 



1 — v 

as X ^ — cxD. It follows that the boundary contribution at x = a 
vanishes also as a — > —00, hence, in the limit as a— > —00, b ^ +00, 

f+00 / - \ f+00 I /|2 

3f?e(A) / [ ^—\u\'^ + \v\'^ ] dx + ^—^dx = 0. (84) 

J -00 \l-v J 7_oo 

But, for ?R.eX > 0, this implies u' = 0, 01 u = constant, which, by 
u{—oo) = 0, implies u = 0. This reduces ()80ap to v' = Xv, yielding 
the explicit solution v = Ce'^^ . By v{-^oo) = 0, therefore, v = for 
3?eA > 0. It follows that there are no nontrivial solutions of (I80p . (I8ip 
for 3f?eA > 0. 

Remark 7. The above energy estimate is equivalent to multiplying the 
system by the special symmetrizer ( i - then taking the 



,0 v/{l 

inner product with (t',u)"^. The analog of the high-frequency es- 
timates of Appendix [XI would be obtained using the alternative sym- 
metrizer Q ^ optimized for its effect on second-order derivative 

term u" /v. This may clarify somewhat the strategy of the energy es- 
timates used in [34112]. 
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C. Quantitative conjugation estimates 

Consider a general first-order system 

W' = A{x,X)W. (85) 

Proposition 5 (Quantitative Gap Lemma |13|,I47| ). Let and 

be an eigenvector and associated eigenvalue of A^(\) and suppose 
that there exist complementary generalized eigenprojections (i.e., A- 
invariant projections) P and Q such that 

< (j^^-f,x ^ > (86) 

\{A-A+){x)\<C2e-''^ x>0, 

with < fj < rj. Then, there exists a solution W = e'^^^V{x, X) of 
(I85D with 

provided {rj — fj)^^CiC2e'~'^^ < e. 

Proof. Writing V' = {A+ — fi'^)V + {A — A-^)V and imposing the 
limiting behavior V{+oo,X) = , we seek a solution in the form 
V = TV, 

I- + 00 

TVix) := y+ - / Pe(A+~M+)(^-?/)(^ _ A+)V{y)dy 

J X 

+ r Qe^^+-^^^^^-y\A-A^)V{y)dy, 



from which the result follows by a straightforward Contraction Map- 
ping argument, using ([86]) to compute that 



/■ + 00 

\TVi - TV2\{x) = - / PeiA+~^^^)i-y){A - A+){yi - V2){y)dy 

J X 

+ £ Qe(^+-'^^)(^-^)(A - A+)(Vi - V2){y)dy 

r+oo 

< C1C2 e-^^^^-y^e-^ydyWVi - V2\\l^[l,+oo) 

— Wl - ^2 L°°fL,+oo)) 

and thus llTl^i - TV2\\lo^[l,+oo) < ^^f^ll^i " ^2||l-[l,+oo)- 



Stability of isentropic viscous shock profiles 



35 



Corollary 5. Let and be an eigenvector and associated eigen- 
value of Aj^{X), where A-^. is n x n with at most k stable eigenvalues 
and 

max\{A+ - < Co; | - < C'se"''^ x > 0, (88) 

< 17 < r/. Then, there exists a solution W = e'^^^V{x, X) of ([85]) 
with 

\V+{X)\ - 6-{v-v){l-e) ^ - ' ^ ^ 
^■=2E^' provided < e. 

Proof. Without loss of generality, take /i = 0. Dividing [—r], —fj] into 
k + 1 equal subintervals, we find by the pigeonhole principle that at 
least one subinterval contains the real part of no eigenvalue of A^. 
Denoting the midpoint of this interval by —fj > fj, we have 

uim\?flea{A+)-fi\>6:=^^. (90) 

Defining P to be the total eigenprojection of A^ associated with 
eigenvalues of real part greater than rj and Q the total eigenprojection 
associated with eigenvalues of real part less than 57, and estimating 
Pe^+^, Qe^+^ using the the inverse Laplace transform representation 



2m 



— / e^^X-z - A+)-^dz, 
m Jr 



with r chosen to be a rectangle of side 4nCo centered about the 
real axis, with one vertical side passing through SReA = —fj and the 
other respectively lying respectively to the right and to the left, and 
estimating 

\{X-A+)-^\ < nlC^-^S-"" 

crudely by Kramer's rule, we obtain (j86|) with Ci = 16nn!CQ (5""", 
whence the result follows by Proposition [5l 

Corollary 6 (Quantitative Conjugation Lemma). Proposition 
[2 holds with 

r.,T N , 1- N 16nn!(Co)"C2e-''^ 

C(L, r?,r/, max I (M+),j I, dim M+) = . 

n := (dimM+)2, k := (<^^-^M^)^-dimM ^ ^^^^ < e. 
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Proof. Writing the homological equation expressing conjugacy of variable- 
and constant-coefficient systems following |35j . we have 

P' = M+P - PM+ + 0M. 

Considering this as an asymptotically constant-coefficient system on 
the n^-dimensional vector space of matrices P, noting that the linear 
operator A^+P := M+P — PM+, as a Sylvester matrix, has at least 
n zero eigenvalues and equal numbers of stable and unstable eigen- 
values, we see that the number of its stable eigenvalues is not more 

2 

than k := " ^" , whence the result follows by Corollary [SI 
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